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Abstract 

We present a microscopic calculation of the ^He /3-decay into the ground state of ^Li. To this 
end we use the impulse approximation to describe the nuclear weak current. The ground state 
wave functions are obtained from the solution of the nuclear 6-body problem. The nucleon- 
nucleon interaction is described via the J-matrix inverse scattering potential (JISP), and the 
nuclear problem is solved using the hyperspherical-harmonics approach. This approach results 
in numerical accuracy of about 2 per mil in the transition matrix element. Bearing in mind that 
the contribution of meson-exchange currents to the transition matrix element is about 5%, these 
results pave the way for accurate estimation of their effect. 
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1. Introduction 

The nuclear weak current contains a well established leading 1-body terms, which at 
low energy are known as the Fermi (F) and Gamow- Teller (GT) operators, and sublead- 
ing two- and many-body currents, generally known as meson-exchange currents (MEG), 
which are a subject of vast contemporary research. The lightest nucleus that undergoes a 
/3-decay is the triton. However, the theory of nuclear weak interaction cannot be checked 
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in the triton since its half-life is used to adjust a free parameter in the axial MEC. There- 
fore, the lightest nucleus that can provide a test to the theory is ^He. ^He (J'^ = 0+) is 
an unstable nucleus, which undergoes a (3 decay with a half-life T1/2 = 806.7 ± 1.5 msec 
to the ground state of '^Li (J'' = 1+) [IJ. 

So far, a microscopic calculation of ^He from its nucleonic degrees of freedom, failed 
to reproduce the /3-decay rate. A comprehensive study, accomplished by Schiavilla and 
Wiringa [2], has used the realistic Argonne wl8 (AVIS) nucleon-nucleon (NN) potential, 
combined with the Urbana-IX (UIX) three- nucleon- force (3NF), to derive the nuclear 
wave functions, through the variational Monte-Carlo approach (VMC). The model used 
for the nuclear weak axial current includes one- and two-body operators. The two-body 
currents are phenomenological, with the strength of the leading two-body term - asso- 
ciated with A-isobar excitation of the nucleon - adjusted to reproduce the GT matrix 
element in triton /3-decay. The calculated half-life of ®He is over-predicts the measured 
one by about 8%. An unexpected result of the calculation, was that two-body currents 
lead to a 1.7% increase in the value of the Gamow- Teller matrix element of ^He, thus 
worsening the comparison with experiment. The authors of this paper have presumed 
that the origin of this discrepancy is either the in approximate character of the VMC 
wave functions, or in the model of the weak nuclear current. 

The /3 decay rate of ^He is proportional to the square of the GT matrix element. The 
difference between the one-body contribution to the ^He-^Li GT matrix clement and the 
experimental value, 2.173 ± 0.002, is of the order of few percent. A result which on the 
one hand is very satisfying, but on the other hand implies that numerical accuracy at a 
per mil level is required to assign the ^He /3-decay as a test case for validating the MEC 
model. 

The current contribution is dedicated to explore exactly that point. Having in mind the 
required level of convergence we use the JISP16 potential [3]. The JISP16 NN potential 
utilizes the J-matrix inverse scattering technique to construct a soft nuclear potential, 
formulated in the harmonic oscillator basis, that by construction reproduces the NN 
phase shifts up to pion threshold and the binding energies of the light nuclei A < 4. 

We use the Hyperspherical-Harmonics (HH) expansion to solve the Schrodinger equa- 
tion. The HH functions constitute a general basis for expanding the wave functions of 
an v4-body system [4|5I . In the HH method, the translational invariant wave- function is 
written as 

* = II C[K]nRnip)y[K] S,, U) (1) 
[K]n 

where p is the hyperradius, and Rn{p) are a complete set of basis functions symmet- 
ric under particle permutations as p^ = ji'^'i ~ '"i)^- The hyperangle, ft, is a set 
of 3 A — 4 angles, and y^j^^ifl, Si,ti) are a complete set of antisymmetric basis func- 
tions in the Hilbert space of spin, isospin and hyperangles. The functions y^x]^^, Si,ti) 
are characterized by a set of quantum numbers [K] |6|7j . and possess definite angu- 
lar momentum, isospin, and parity quantum numbers. They are eigenfunctions of the 
hyperspherical, or generalized, angular momentum operator K^, K'^y^K]{^, Si,ti) = 
KiK + 3A-5)y[K]i^,s,,U). 

Our results for the ground state properties of ^He and ^Li, are presented in table [T] 
The energies and rms matter radii are given as a function of Kmax ~ the limiting value 
of K in the HH expansion. As we have used the bare interaction in our calculations 
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the results are variational. Using the formula E{Kmax) — E^o + A/K"^^^ for Kmax > 8 
we have extrapolated the binding energies to the limit Kmax — > oo. It can be seen 
that the extrapolated binding energies are rather close to the experimental values. In 
the last column of the table we present our ^He-^Li GT transition matrix element in 
the impulse approximation, i.e. at the 1-body level. It can be seen that the convergence 
pattern of the matrix element is not regular. Extrapolating its value using the expression 
GT(K,nax) = GToo + B/K^,, for if,„,, > 6, we get GT^o = 2.228(3). This result is 
in accordance with the values GT = 2.28 for AV87TM'(99) and GT = 2.30 for AV8' 
obtained by Navratil and Ormand 0, GT = 2.28 for the N^LO NN-force of Navratil and 
Caurier [9], GT ^ 2.25 for AV18/UIX of Schiavilla and Wiringa [2], and GT = 2.16-2.21 
for AV18/IL2 by Pervin et al. [10]. Moreover, it can be seen that our accuracy estimating 
the GT matrix element is in the level of few per mil. Such an accuracy enables us to 
disentangle numerics from physics and utilizes the ^He /3-decay as a test ground for axial 
MEG models. 

Table 1 

The ^He, ®Li binding energies, rms matter radius, and GT matrix element in the impulse approximation. 
Results for the JISP16 NN interaction. 
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2.165 






Shirokov ei al. |3] 


28.32(28) 


31.00(31) 








Exp. 
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